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Abstract. 

The early-R stars are carbon-rich K-type giants. They are enhanced in ^^C, ^"^C and ^''N, have approximately 
solar oxygen, magnesium isotopes, s-process and iron abundances, have the luminosity of core-helium burning 
stars, are not rapid rotators, are members of the Galactic thick disk and, most peculiarly of all, are all single stars. 
Conventional single-star stellar evolutionary models cannot explain such stars, but mergers in binary systems have 
been proposed to explain their origin. 

We have synthesized binary star populations to calculate the number of merged stars with helium cores which 
could be early-R stars. We find many possible evolutionary channels. The most common of which is the merger of 
a helium white dwarf with a hydrogen-burning red giant branch star during a common envelope phase followed by 
a helium flash in a rotating core which mixes carbon to the surface. All the channels together give ten times more 
early-R stars than we require to match recent Hipparcos observations - we discuss which channels are likely to 
be the true early-R stars and which are not. For the first time we have constructed a viable model of the early-R 
stars with which we can test some of our ideas regarding common envelope evolution in giants, stellar mergers, 
rotation, the helium fiash and the origin of the early-R stars. 



1. Introduction and History 

The early-R stars are one of several classes of carbon 
star, all distinguished spectroscopically by the presence 
of strong absorption due to an excess of molecular car- 
bon. While many of these carbon stars are now relatively 
well understood, the early-R stars have defied satisfactory 
explanation for nearly one hundred years. In a nutshell, 
the problem is to explain how a single star which is a gi- 
ant, and hence has finished core-hydrogen burning, but is 
not sufficiently luminous to have completed core-helium 
burning, can have a surface which is contaminated by ex- 
cess carbon. It may be argued that such a straightforward 
statement of the problem is too simplistic; it probably is. 
Therefore it is necessary to review here the main charac- 
teristics of the early-R stars, to define the problem in more 
detail, and thence introduce one approach to its solution. 



1.1. Spectral class definition 

Stars wit h pecu l iar ba nding in their spectra were iden- 
tified by ISecchil dlSGSl) and were later identified as car- 
bon stars by Rufuj ( 19 id ). They are easily identified by 



strong absorption features due to C2, CN or CH molecules. 
A subset of these, the R-typ e stars, was first classified 
by Fleming fc Pickerind (Il908h based on the observations 
of Fleming fc PickeringI (1890). The R stars "contain rays 
of much shorter wave length than ordinary fourth type 
stars" which, in modern parlance, means they are bluer 
and hotter (Teff ^, 4000 K, similar to the K-type stars, 
e.g. lOominvl . 1 19841) than the normal N-type carbon stars 
(with Toff ^ 3500 K) whi ch are probably asymp totic gi- 
ant branch (AGB) stars (jlben fc Renzinil . 119831 ). The R 



stars are very com mon, apparent ly accounting for 1% of 
all K and G giants ( Dominyl . 19841 ). and are 10 times more 



numerous than N stars according to iBlancol ()l965l l. 
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Later, as part of th e Henry Draper catalogue 
( Camion fc PickeringI 1918f ) and subject to revision by 
Shand ( 1928f ). the R class was split into RO to R8, where 
RO-4 (the hot/early/warm-R stars) are warm, equiva- 
lent to normal K-type stars, while the cool/late-R stars 
of type R5-8 are the carbon-rich equivalent of M stars. 
All R and N stars were rn erged into a single type C 
by iKeenan fc Morgan ( 1941), b u t this was a "retrogres- 
sive step" according to lEggen (|l972l ). who reclassified 
late-R stars as N stars and left the early-R stars as 
a distinct class. Other types of carbon stars compli- 
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cate the issue, especially in binaries where mass trans- 
fer can pollute the secondary star with enough carbon 
to turn it into a dwarf carbo n star or CH star (see 
e.g. Wallerstein fc Knappl . [l995 sections 3 and 5). The J 
stars, which are similar to N stars with enhanced ^'^C, are 
also of uncertain origin but are cle arly redder and mor e 
luminous than the early-R stars ijAbia fc Isernl . l2Q00l ). 
Example spectra of all the carbon star types can b e 
found in the atlas of iBarnbaum. Stone fc Keenan ( 1996h . 
There is an alternative cla s sificati on scheme developed by 
Bergeat. Knapik fc Rutilvl (|2002a[ ) who identify RO-3 stars 
as members of their HC1-HC3 class (where HC stands for 
hot carbon). 

Standard single-star evolutionary models of 
AGB stars undergoin g thi rd dredge-up (e.g. 



by mass transfer to a main sequence companion from an 
AGB star which has undergone third dredge-up. They are 
rich in s-process elements and carbon which were formed 
in the primary AGB star and are all binaries. The pol- 
luted secondary evolves to the C HeB phase, perhaps re - 



maining carbon- rich (although see IStancliffe et al.1 . l20Q7n 



Metallicity affects the CH-formation process in two ways. 
First, third dredge-up is m ore efficient at sub-so lar metal- 



licities compared to solar i Karakas et al. . 2002f ) 



Karakas. Lattanzio fc Pols . 2002[ ) correctly, if 



some- 



primary carbon is made and transferred to the main se- 
quence star. Second, the initial abundance of oxygen is 
lower than at solar metaUicity so less carbon is required 
to be mixed into the stellar envelope to form a carbon 
star (for which C/O > 1 by number). There is a threshold 
metallicit y, around Z ^ OAZ p, , above which CH stars can- 
not form I Abia et "all . 1200^ . Both these facts are at odds 
with the properties of the early-R stars, which have solar 
or sli ghtly sub-solar iron and oxygen abundances ijPominyl . 
19841 ). 



1.4. Binary Fraction 



times only qualitatively, predict most of the properties 
of the late-R and N type carbon stars. The early-R stars 
remain a distinct group of stars with an unknown origin, 
just as they were in 1908. More information abou t the 
carbon star family is given by I Abia et all ll2003h. a l- 
though also useful ar e the reviews ofjUovd Evanf^ ('1986'), 

jodders fc FeglevI ( 1998h and IWallerstein fc Knapp The nai l in the AGB rn ass-transfer cofiin came with the 
IQQSl l. work of lMcClurd l|l997t ) who found that 22 early-R stars 

- all those in his sample - are single stars, even though 
20% of late-type giants are binaries. The implication is 
that the early-R stars originate in binary mergers, as it is 
hard to envisage a physical process which makes carbon 
stars only in single stars and not in wide binaries. That 
they are single and not s-process enhanced h as been used 



1.2. Spectroscopic Studies 

Spectroscopic studies of the early-R stars have shed some 
light o n their or i gin. T he most comprehensive analysis is 
that of Dominvl (1984). He found that early-R stars, de- 
spite having C/O > 1, have a low ^^C/^^C ratio (sim- 
ilar to the J stars), a solar or slightly sub-solar iron 
abundance, enhanced nitrogen relative to solar, a so- 
lar oxygen abundance, solar ^^'^*0/^^0 ratios and no s- 
process enhancement. This is in contrast to N-type car- 
bon stars, and late-R stars, which show evidence of s- 
process elements, supporting the i dea that late-R stars 
are really just misclassified N stars. Zamora et al. 1 20061 ) 
confirm the findings of Dominy and support the idea 
that late-R and N type stars are probably equivalent. 
Late-R stars pulsate, like N stars, but early-R sta.rs do 
not, like K-type stars dBarnbaum, Stone fc Keenaiil . 1996t 



to di s tinguish e arly-R stars from CH stars IjAbia et al 



|2003[ ) . iMcClurj also found that the early-R stars are not 
rapidly rotating, a feature which must be explained by any 
potential early-R star model. 

1.5. Space Density and Distribution 



Bergeat. Knapik fc Rutilvl . l2002bl fig. 8). 

The luminosity of early-R stars is around 100 
typical of core helium burning (CHe B) stars (the red 
cluni p in a col o ur-magnitude diagram, Cannon fc Llovi 
1969[ ICannonl . 1 1970l : (Faulkner fc Cannonl . Il973l : IScaL 
19761 ) ■ while the luminosity of N and l ate-R stars i s more 
than 1000L(7i, typica l of AGB stars d Vandervort . 1958t 
Baumertl . Il974t IScalol . Il976t iBergeat et all . l2002al l. 



Knapp. Pourbaix fc JorissenI (|20011 ) determined the abso- 
lute magnitudes and space density of the early-R stars 
using Hipparcos data. They found that the early-R stars 
have magnitudes similar to the red clump (i.e. CHeB) 
stars, the space density of the early-R stars is 4.5 — 
15 X 10~^pc~ ^ and the R to red cl ump number ratio is 
0.04 - 0.14%. lBergeat et all l|2002b[ ) calculated the space 



1.3. CH Stars 

There remains the possibility that the early-R stars are 
the CHeB counterparts of the CH stars in which, for 
some reason, s-process elements are undetectable or not 
present. CH stars, both dwarfs and giants, are formed 



density in the Galactic plane to be 1.66 x lO^^pc"^ for 
early-R star^, 13.4 x 10~®pc~^ for N stars, a ratio of 
~ 16. The two papers give different number densities, but 
actually their projected number densities are about the 
same. 

That t he R st a rs are Galactic d i sk obj ects was recog- 
nised by lEggenl (|l972f ). IPominvl (|l984h finds that N 
stars are younger disk o b jects t han, and distinct from, 
the R stars, while Ishida Jl960l) . StephensonI 1 19731 ) and 
Barbaro fc Dallaportal ()l974l ) found that the N and R 



^ Strictly, iBergeat et aP l|2002bl ) define new spectroscopic 
types, so we assume, as they state, that early-R stars are their 
HC type, N stars their CV type. 
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stars are distributed differently across the sky, with N stars 
more condens ed in the Galac t ic plan e. These results were 
confirmed by iBergeat et al.1 I 2002bl ) who found that R 
stars are, on average, three times further from the Galactic 
plane. The solar or slightly sub-solar iron abundance and 
velocity dispersion of the R st ars suggests they are mem- 
bers o f the Galactic thick d i sk llSanfordl. 1944 1 Vandervort , 



Recently, full three dimensional models 
of the helium flash have been constructed 
Dearborn. Lattanzio fc Eggletonl . 20061 : Lattanzio et al 



Bergeat et al 



2nn2bh " 



1958t iDahnl. Il964l: lEggenl . Il972t iKeenan fc Heck 



19941 



Third dredge up only occurs in 
metallicity AGB stars above abou t 
1.5 M^^ llBessell. Wood fc Llovd Evansl . 



solar- 
1. 3 - 
1983t 



although these are of single stars, not binary merg- 
ers (in the latter paper slow rotation was introduced). 
Simulations of HeWD mergers have been constructed 
with a smooth particle hydrodynamic s (SPH) approach 
( Guerrero. Garcia-Berro fc Isernl . 2004^ . While these do 
not evolve to the helium flash, they do show that the 
cores lose very little mass, or angular momentum, during 
their merger. 



Karakas. Lattanzio fc Pols[ 2002f ). corresponding to J. 7. Summary 



a stellar lifetime of 5 Gyr. Hence the early-R stars cannot 
be intrinsic AGB stars if they are older than about 5 Gyr. 



1.6. Stellar Models 

The early-R stars pose a problem for stellar evolution the- 
ory. Standard models of single stars are not carbon rich 
except for thermally pulsing AGB stars and some mas- 
sive Wolf-Rayet stars. R-stars are too dim to be either 
of these, indeed their luminosities are those of the CHeB 
(red clump) evolutionary phase. The problem we face is 
twofold: when is carbon made in these stars and how do 
we get it to the surface? 

The binary star merger model satisfles most of the 
observational constraints. We know there are no early- 
R stars dimmer than the red clump, so the merger must 
lead to a carbon-rich, CHeB star. This naturally implies 
that something odd happens immediately prior to the 
CHeB phase during helium ignition at the tip of the giant 
branch. Canonical models of the helium flash do not pre- 
dict mixing of carbon-rich material from the core to the 
stellar surface ( Harm fc SchwarzschildL 19661 ). However, 
these models are one-dimensional and non-rotating, while 
a binary merger leads to a three-dimensional, rotating 
system, so it is quite possible that the merger and/or 
ignition leads to non-standard mixing. One-dimensional 
models with parameterised rotational mixing and/or low 
metallicity (with off-centre helium ignition) have led to 
some mix ing of carbon-rich core r naterial with the stellar 
envelope jMengel fc Gross j_ 1976 : Paczynski fc Tremaind . 



19771: iHollowell. Iben fc Fuiimotol . Il990l l. The models of 
IbenI (|l990l ) simulate accretion of helium onto a helium 



white dwarf (HeWD) and show that off-centre helium ig- 
nition is a natural consequence of this process. 



A 
been 



number 
carried 



of two-dimensional simulations have 
out by Deupre e and collaborators 



JCole fc Deupred. Il980l . Il98ll : lOeupree fc Wallacd 



Il987t Deupree . 199w some of which suggest the possi 



bility of mixing material from a core helium flash into 
the stellar envelope. It is hard to draw a conclusion from 
these models as the results vary with both resolution and 
model sophistication, and it is not clear how to apply the 
results to a rapidly rotating binary merger. 



In this paper we will make, for the first time, a quantitative 
estimate of the number of and properties of early-R stars 
using our binary population nucleosynthesis model. We in- 
vestigate possible channels for early-R star formation by 
a binary merger process and subsequent heHum fiash in 
a rapidly rotating core and also extrinsic (CH-star) chan- 
nels by accretion from a companion. We determine the 
effect of varying model parameters on the early-R to red 
clump and early-R to N star ratios, and compare these 
to observations. Section [2] describes our model, section [3] 
our results, section [3] discusses the ramifications and we 
conclude with some ideas for future research. 

2. Modelling the Mystery 

We model populations of si ngle and binary stell a r star s 
with the synthetic model of lHurlev. Tout fc Polsl ll2002l) 
updated to i nclude nucleosynthesis ijlzzard et al.l . 12004 : 
Izzard et al . 20061 and changes outHned below). Our 
model follows the evolution of stars using analytic fits 
to luminosity, radius, core mass etc. Interactions with 
a companion due to tides, wind accretion and Roche- 
lobe overfiow (RLOF) are taken into account. Common 
envelope evolution is treated with a n a-prescription 
()Livio fc SokeA \l98A iTout et al.1 . Il997[ ). where a is the 
fraction of the orbital energy which is transferred to 
the envelope during the spiral-in phase. In most of our 
models we use a = 1 and an env e lope binding en- 
ergy factor A = 0.5 i Regos fc Tout . 1995l l. Our syn- 
thetic AGB model is based on the full evolutionary mod- 
els of Karakas. Lattanzio fc Polj 1 2002h and includes sur- 
face abundance changes due to third dredge-up and hot- 
bottom burning. We also follow the surface abundances o f 
massive stars based on the models of Drav fc Tout I 2003h . 
as well as yields due to supernovae and novae, but these 
are not important with regard to the early-R stars. 



2.1. Stellar evolution 

Our synthet ic stars ha ve solar-scaled initial abundances 
according to I Anders fc Grevesse (19891 ). with metallicities 
10^'' < Z < 0.03. There are many oth er parameters in ou r 
model which are discussed in detail in Izzard et al. 1 2006[ ). 
Where a parameter choice affects our results, we vary the 
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parameter within a reasonable range (see section lS^ . The 
main parameters which influence R-star formation are the 
metallicity Z and the common envelope parameters a and 
A. These are Z = 0.02 (which we refer to as solar metal- 
licity), a = 1 and A — 0.5 unless stated otherw ise. 

We have updated the llzzard et al. ( 20Q6f ) treatment 
of case-B RLOF in binaries with primary masses be- 
tween 0.75 and 2.0 M©, which turn out to be our pro- 
totype early-R stars. We follow the surface abundances 
of ^H, '*He, ^^C, ^''N and ^^O as a function of mass co- 
ordinate based on detailed terminal main-sequence mod- 
els constructed with the TWIN st ellar evolution code 
l Eggleton fc Kiseleva-EggletonL 2002l l. As mass is stripped 
during RLOF, transferred material is enhanced in helium 
and nitrogen but deficient in carbon and oxygen because 
of CNO cycling during the main sequence evolution of the 
primary. This is transferred to the secondary star, enhanc- 
ing its abundance of helium and nitrogen while reducing 
hydrogen and carbon. 

We mod i fy the common envelope prescription of 
Hurlev et al. ( 2002f ) by removing their algorithm which 
instantaneously ignites merging degenerate helium core^. 
Instead, in our model the less massive core is dis- 
rupted into a disc around the more massive core and 
accretes until a single , rapidly rotating core is formed 



I Guerrero et al.l . 120041 1. The merged core then grows by 



hydrogen-shell burning until helium ignites in a nuclear 
runaway (the helium flash) . For all our merged stars we as- 
sume that the ignition process, in a rapidly rotating core, 
mixes some carbon int o the stellar envelope as found by 
Paczvnski fc Tremaind (|l977l ). The star then settles into 

its CHeB phase as a n early-R star. 

According to the iHurlev. Pols fc ToutI ()2Q0Q| 1 formal- 
ism, helium ignites when the stellar luminosity reaches 
a critical value which is a function of the stellar mass 
and metallicity. This may not apply to our merged stars 
which have abnormally large cores for their total mass 
and stage of evolution, but given that we are ignoring 
the lifting of degeneracy due to the merger we can prob- 
ably do no better without detailed models of early-R 
star progenitors. In ou r main merge r chan nels (see be- 
low) application of the Hurlev et al. ( 2000f ) prescription 
leads to minimum helium-ignition core masses of between 
0.35 and 0.45 M©, compatible with the off-centre He WD 
ignition models of Ibe which have a core mass of 

0.38 Mq. In contrast, an SPH simulation of the merger of 
two 0.4 Mq HeWDs does not lead to a helium flash even 
though the maximum temperature rea ched is 2 x 10* K 
I Guerrero. Garcia-Berro fc Iserrj . 2004 ) - any excess en- 
ergy goes into lifting the degeneracy and expansion of the 
white dwarfs. These simulations do not model the evolu- 
tion of the star beyond a few minutes after the collision, 
so it is possible that the core contracts and/or increases in 
mass on a longer timescale, but ignites while still rapidly 
rotating. 



2.2. Stellar population model 

Our single star and binary primary masses are dis- 
tri buted according to the ir iitial mass function (IMF) 
of iKroupa. Tout fc Gilmord ()l993h in the rangcH 0.1 < 



M/Mq < 80, secondary star masses are chosen from 
a distribution which is flat in g = M1/M2 such that 
O.IMq/A/i < <7 < 1 and initial separations a are cho- 
sen from a distribution flat in log a for 3 < a/R© < 10^. 
We usually assume all binary orbits are circular (e = 0) 
and to a maximum age of 13.7 Gyr. Given these distri- 
butions, the probability assigned to each binary system is 
^/{Mi, M2, a) and the contribution to the number of stars 
of a given type is /^;^^ lM^,M2,a ^ ^ '^(type) x * dt, where 
(5(type) = 1 when the star is of the required type, oth- 
erwise. A similar calculation is performed for single stars 
and the results combined with a 50% binary fraction. We 
set S, the star formation rate, to 1 because we compare 
only number ratios and relative number counts. 



2.3. Age selection criterion 

The early-R stars are associat ed with the Galacti c 
disk, particularly the thick disk i Bergeat et al. . 2002bh . 
We select model stars older than the lower 



limit 

of the age-metallicity relation for t hick disk stars of 
Bensbv. Feltzing fc LundstromI ()2004l l. which we flt to 
5.61 - 6.68/ Gyr where / = max (-0.65, [Fe/H]). At so- 
lar metallicity ([Fe/H] = 0) this gives a minimum age of 
5.61 Gyr, which is too old for intrinsic carbon star forma- 
tion according to our models i Karakas et al. . 2002f ). 



The imposition of the age limit is equivalent to ending 
star formation 4.5 Gyr ago, which is not applicable to the 
whole Galactic disk, even if it is correct for the thick disk. 
It implies, for example, that there are no intrinsic N type 
stars, which is incorrect in the thin disk so we must be 
careful when comparing R to N star number count ratios 
with the results of our model. 



2.4. N and red clump stars 

We deflne N-type carbon stars as giant branch (GB) or 
AGB stars with Nq/Nq > 1. All our N stars are extrinsic: 
there are no intrinsic carbon stars in our simulated pop- 
ulations because such stars are younger than our ~ 4Gyr 
age limit. 

The red clump is synonym ous with the CHeB p hase of 
stellar evolution (as deflned by Hurlev et al. , 2002f ) in low- 
mass stars, equivalent to the horizontal branch at higher 
mass. 



^ The energy released in their algorithm usually disrupts the 
star, which we consider unrealistic. 



^ After some low-resolution trial runs the primary mass range 
was reduced to 0.3 < M/Mq < 8 with little change to our 
results. 
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2.5. Early-R star formation channels 

The iHurlev et all ll2nn2h model defines many common- 
envelope merger channels which lead to a CHeB star, any 
of which could be the early-R stars. We define our R-star 
formation channels as follows: 

Rl. Merger of a naked main-sequence helium star 

(HeMS) with a GB star. 
R2. Merger of a core-helium burning star (CHeB) with 

a GB star. 

R3. Merger of a helium white dwarf (He WD) with the 

core of a GB star. 
R4. Accretion from an AGB companion, which is now 

a white dwarf (at low metallicity, these are CH stars). 
R5. Merger of a HeWD with a Hertzsprung gap (HG) 

star which has a partially degenerate helium core. 
R6. Merger of a GB star with a HG star. 
R7. Merger of two HG stars. 

R8. Merger of an AGB star with a helium-cored star. 
The helium and CO cores mix completely to give a 
new, CO-rich helium core which behaves as an evolved 
CHeB star. 
R9. Merger of two naked helium stars. 
RIO. Merger of a CO or ONe white dwarf (i.e. the core 
of an AGB star) with an AGB star. This is similar to 
channel 8 but one star has its envelope stripped. 
Rll. Merger of two GB stars, similar to channel 3 but 
both stars have a hydrogen envelope. 
Our early-R stars are defined as those which have gone 
through one of the above channels, are of the appropriate 
age (see section 12. Sp and are in the CHeB (red clump) 
phase. We assume that all the above mergers make early- 
R stars by the rotating core helium fiash mechanism. 

3. Results 

Table [T] gives the relative formation rate and numbers 
of stars which form in each of our model early-R star 
channels in our standard Z = 0.02 population. The R3 
(HeWD-GB) and R6 (HG-GB) channels together repre- 
sent more than 90% of our early-R star progenitors - we 
discuss these stars below in some detail. We identify chan- 
nel R4 with the CH stars. These form in appreciable num- 
bers at low met allicity, but not a t all at solar metallicity, 
as predicted by Abia et aU ()2002l ) - we discuss these sep- 
arately in section [331 R5 (HeWD-HG) and Rll (GB-GB 
merger) possibly are early-R stars, forming most of the 
remaining 10%. There are small contributions from the 
R7 and R8 channels and none in the Rl, R2, R9 and RIO 
channels. 

The parameter space of initial masses and periods 
which goes on to form our early-R star candidates is shown 
in figure [TJ 



3.1. Subtype R3: HeWD + GB merger (Z = 0.02j 

The R3 channel consists of HeWD-GB mergers and dom- 
inates our early-R star progenitor population at solar 



Channel 


Formation Rate 


Number Fraction 


R3 


0.766 


0.705 


R5 


0.016 


0.041 


R6 


0.169 


0.196 


R7 


0.002 


0.005 


R8 


0.014 


0.004 


Rll 


0.033 


0.048 



Table 1 . The fraction of each merger channel which con- 
tributes to our synthetic, standard Z = 0.02 early-R 
star population as defined in section [231 The left column 
(Rl. . . 11), the middle column the relative formation rate 
for each channel and the right column the relative number 
of early- R stars (assuming a. constant star formation rate 
and the Kroupa et al . 19931 IMF). The missing channels 
do not to our R star population and channel R4 is not 
a merger (we identify these with the CH stars and treat 
them separately, see section [3731 . 



metaUicity (Z ~ 0.02). We have split the R3 type into sub- 
types R3a and R3b representing initially short and long 
period binaries as follows: 

- R3a(77%ofR3): Afi = l-1.7Mo,M2 = 0.3-1.0 Mq, 
P = 1 - 2 days (5 < a/R© < 8) , 

- R3b (23% of R3): Afi = 1 - 1.7 Mq, Ma = 0.95 - 
1.3 Mq, P = 150- 700 days (150 < o/R© < 450) . 

Both the R3 subtypes originate from two phases of mass 
transfer, the second of which is the common envelope 
phase which leads to merger of two helium cores. 



3.1.1. R3a Sub-type 

In our R3a stars the first phase of mass transfer is conser- 
vative RLOF when the more massive primary exhausts its 
core hydrogen, crosses the Hertzsprung Gap (HG) and be- 
gins ascent of the giant branch (GB). Initially, while <? < 1, 
the orbit shrinks, but when enough mass is transferred 
that q> I the orbit widens. The primary completely loses 
its envelope and is l eft as a ^ 0.2 M© HeWD (similar to th e 
detailed models of Kippenhahn. Kohl &: Weigertl . Il967 ). 
The secondary accretes a significant amount of mass and 
becomes a blue straggler because it is more massive, hence 
bluer, than it should be for its age. It is also nitrogen and 
^^C rich because it accretes material which is stripped 
from CN-processed layers deep inside the primary. The 
separation of these binaries, prior to the evolution of the 
secondary up the giant branch, is about 20 R©. 

The R3a channel contains the majority, 77%, of the R3 
channel stars because they have a higher formation rate 
and they have smaller cores than R3b so spend longer in 
the CHeB phase. 



3.1.2. R3b Sub-type 

The R3b stars begin their evolution in a relatively wide bi- 
nary. When the primary evolves onto the GB it overfiows 
its Roche lobe and a common envelope results. During 
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Fig. 1. (a) Initial secondary mass vs initial primary mass and (b) initial period vs initial primary mass for all our 
early-R star merger progenitors (i.e. not including channel R4, the CH stars). {Z = 0.02, a = 1, A = 0.5.) The symbols 
are as follows: R3 (HeWD-GB) small grey circles, R5 (HeWD-HG) squares, R6 (GB-HG) upward-pointing triangles 
(A), R7 (HG-HG) downward-pointing triangles (v), R8 (AGB) large diamonds (o), RU (GB-GB) open circles. The 
contribution from omitted channels is negligible. In panel (b) it is clear that there are two distinct populations of R3, 
R6 and RU, the short period binaries which have one common envelope phase, and the long period binaries which 
have two. 
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Fig. 2. Pre-common envelope HeWD mass vs giant mass (a) and giant core mass (b) for our R3 channel. {Z = 0.02, 
a = 1, A = 0.5.) In both plots the group of points on the left side comes from our R3a (initially short-period binary) 
models, while the points on the right are from the (initially long-period) R3b channel. The R3b stars have such large 
cores that they are probably too bright to be the observed R stars. 



the spiral-in of the heHum core and main sequence star, 
the envelope of the giant is lost, leaving a 0.38 — 0.48 Mq 
HeWD with a 1 — 1.2 Mq main sequence (MS) star in a 
reasonably close (a ^ 20 R©) binary. The HeWD is about 
twice as massive in this scenario as compared to R3a be- 
cause the giant evolves further up the giant branch before 
the first mass transfer. Note that in this scenario, the sec- 
ondary accretes only a small amount of mass (typically 
~ 0.01 M0) so while technically it is a blue straggler, it 
might not be detected as such. Also, the secondary does 
not accrete much ^'^C or ^''N from the stripped primary. 

3.1.3. Formation of R3 stars 

In both the R3a and R3b channels, after the first mass 
transfer the secondary evolves on its nuclear timescale and 
eventually exhausts its core hydrogen. As it crosses the HG 
and ascends the GB, RLOF begins and common envelope 
evolution results. 

Figure [2] shows the distribution of masses and core 
masses just prior to the final common-envelope phase and 
core merger. In both panels of the figure the R3a stars 
are on the left and the R3b on the right. The giant core 
and the HeWD merge during the ensuing common enve- 
lope phase to form a new single GB star with a core of 
0.3 - 0.45 Mq for the R3a and 0.5 - 0.7 Mq for the R3b 
channel. The post-merger object is not immediately an 
early-R star, but is a GB star with an abnormally large 
core. In the case of R3b and more massive R3a stars, the 
core is massive enough that the star very quickly ignites 
helium. In the lowest-mass R3a stars some hydrogen shell 
burning drives the core mass up to a minimum of 0.36 Mq 
when heHum ignites. Once helium ignition has started, 
the luminosity of the star drops and it settles into the red 
clump as a candidate early-R star. 



Because the R3b stars have massive cores they have a 
mean luminosity of about 650 Lq. This is rather large, 
both compared to the observed early-R stars and to 
our R3a channel, which has a mean luminosity of about 
170 Lq. Once the core-heHum burning R-star phase is 
complete, these stars ascend the AGB, possibly as J-type 
stars (see section [421) • 

3.2. The other merger channels (Z = 0.02) 

The other merger channels arise in much the same way as 
R3, with initially close binaries undergoing conservative 
mass transfer onto a companion and initially wide bina- 
ries passing through two common-envelope phases. These 
other phases are considerably rarer because the initial bi- 
nary parameters must be just right such that the second 
mass-transfer phase occurs e.g. during the HG rather than 
GB phase. Consequently, channels R5 (HeWD-HG), R7 
(HG-HG), R8 (AGB) and Rll (GB-GB) are rare com- 
pared to R3. The R8 (AGB-merger) channel is almost cer- 
tainly not related to the early-R stars because the result- 
ing merged stars have luminosities around IO'^Lq, typical 
of stars with evolved CO cores (late-R stars or N stars). 
We make no R9 or RlO stars at solar metallicity and al- 
though some RlO (CO WD-GB or ONe WD-GB mergers) 
systems do exist at lower metaUicity their numbers are 
very few. We never make any Rl or R2 systems because 
we select only old stars and potential R2 stars merge as 
GB-GB systems (i.e. prior to the GB-CHeB phase). 

Channel R6, the merger of a GB star with a HG star, is 
relatively common. Most of these systems (85%) are sim- 
ilar to the R3a systems but with an initially more mas- 
sive secondary, such that it is a HG star (rather than a 
MS star) when the primary overfiows its Roche lobe. As 
a result, a common envelope forms during the first mass 
transfer phase and the cores merge. 
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Fig. 3. Our model predictions for the relative fractions of 
the R3 and R4 channels as a function of metallicity. While 
R4 dominates at low metallicity (these are halo CH stars), 
channel R3 (HeWD-GB mergers) dominates at around so- 
lar metallicity - these are our prime candidates for the 
early-R stars. 




R3 (HeWD GB) R6 (GB HG) A R8 (AGB) — 0- 

R5 (HeWD HG) □ R7 (HG HG) V Rl 1 (GB GB) —Q- 

Fig. 4. The relative number of our early-R star merger 
candidates as a function of metallicity, arbitrarily normal- 
ized to 1 for the R3 channel at Z = 0.03. The frequency 
of most merger types decreases as the metallicity drops. 
The symbols are the same as in figure [H 



The remainder of the R6 systems, and most of the Rll 
systems, are initially wide binaries containing a HeWD 
and a HG (or GB, in the case of Rll) star. The wind from 
the HG star is accreted onto the HeWD at a rate sufficient 
to form a new envelope and rejuvenate the HeWD as a 
GB star. A common envelope phase soon follows and the 
HeWD merges with the HG star. These systems suffer the 
same problem as the R3b channel, their cores are massive 
(0.6- 0.8 M0) and they are are too luminous (^ 600 Lq) 
to be early-R stars. 

3.3. Channel R3 vs R4 (CH) as a Function of 
Metallicity 

Figure[3] shows the relative number of the R3 and R4 chan- 
nels as a function of metallicity. The CH-star channel R4 
dominates our early-R star population for Z < 0.004 while 
above this metallicity the binary merger channel R3 is pre- 
dominant. The reason for this is a combination of increas- 
ing oxygen abundance (proportional to Z) and decreasing 
third dredge-up efficiency, as discussed in section fOl Our 
conclusion is that for Z < 0.004, e.g. in the Galactic halo, 
most carbon-rich red clump stars are probably CH stars, 
with enhanced s-process abundances and a binary com- 
panion, while at higher metallicity (in the Galactic disk) 
the merger model dominates and all carbon-rich red clump 
stars should be single early-R stars. This is just as is ob- 
served (see section [T2)) . 



3.4. The Effect of Metallicity on Merger Channels 

Figure |4] shows the effect of metallicity on the relative 
number of each merger channel. Lowering the metallicity 
reduces the number of mergers of all types mainly be- 



cause of our age-selection criterion (see section 12. 3p . The 
R5 channel (HeWD-HG) increases because as the metallic- 
ity drops the stellar lifetime for a given mass decreases, so 
lower mass stars evolve off the main sequence faster than 
they would at solar metallicity. The IMF favours low-mass 
stars, hence the increase in R5. 

3.5. Early-R to red clump ratio 

The ratio of the number of R to red clump s t ars w as deter- 
mined to be 0.04 - 0.14% bv iKnapp et all l|200ll l. Figure 
[5] compares the results of our models to this ratio and it is 
obvious that our model overestimates the number of early- 
R stars (or underestimates the number of CHeB stars). In 
our R-star selection criteria we have tried to be as inclu- 
sive as possible, so (for example) if only 10% of our R3 
stars go on to mix carbon into their envelopes when they 
ignite helium, then our models would agree with the obser- 
vations - we discuss this further in section HI It is unlikely 
that we have underestimated the number of CHeB stars 
by a factor of ten. 

3.6. The Early-R to N ratio 

The ratio of the number of R to the number of N stars 



is about 10 (jDominvl . 119841 ) . Our models give TVr/A^n 
10 — 100, but this is spurious because we do not include 
young, intrinsic N stars and some of the R stars used to 
calculate the observed ratio are probably cool/late-H (i.e. 
N type). Also, the obs erved rat i o depe nds very strongly on 
Galactic latitude (e.g. Sanfordl . 1944 . com pare Figs 1 and 



2) and /or height above the Galactic plane ijBergeat et al 



2002bl ). To test this one would build a Galactic population 
model, with separate thick and thin disk populations, but 
this is beyond the scope of this paper. 
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Fig. 5. Our calculated R to CHeB (red clump) number ra- 
tio as a function of metallicity for channels R3 to R7. The 
other channels are negligible. The symbols are the same 
as in figure m The thick err o r bar shows the observational 
bounds from iKnapp et al.l ( 2001 ) for the Galactic disk. 
The R3, 5, 6 and 7 subtypes could all be the real early-R 
stars, although in the case of R3 and R6 only a subset 
of our model stars can be true early-R stars because we 
make too many. 



On the other hand, our model results are not in- 
consistent with the observations. We can estimate the 
ratio N^/NcHcB from the observed N^t./NN by c alcu- 



lating Nn/Nf^s x A^n/A^choB- iBergeat et all l|2002bh give 
Nr/N^ ^ 0.32 in the Galactic plan^, while we can ap- 
proximate N^/NcHcB from AiN/A^cHeB ^ 2/150 = 0.013 
where the Ats are the lifetimes (in Myr) of the AGB and 
CHeB evolutionary phases for aZ — 0.02, 1.3 Mq star. 
This results in NTi /Ncyt^,n ^ 0.4% , similar to the 0.1% 
Knapp. Pourbaix fc Jorissera 1 200ll l find and within the 
range of our results. 

3.7. Early-R/ GK giant ratio 



The R to GK giant ratio is 1% according to iDominv 
1 1984 . although the source of this number is not given 
(it is pr esumably ca l culate d from the carbon star cata- 
logue of IStephenso 3, Il973l and probably includes late-R 
stars) . If we define GK giants as GB and AGB stars with 
3800 < Toff/K < 5850 which satisfy our age criterion, 
which is quite conservative as we should probably include 
CHeB stars as well, then at solar metallicity {Z — 0.02) 
our early-R to GK giant ratio is about 0.5%. Given that 
our early-R star count is as high as we can possibly make 
it (it should be a factor of ten less to match the R to red 
clump ratio) and we do not include younger GK giants, 
ou r R to GK giant ratio is clearly different to that quoted 
by IPominyl . 



■* Again, we associate their type CV stars with N stars, and 
their HC stars with early-R stars. 
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Fig. 6. Distribution of the mass of carbon which must 
be dredged up in order to convert our merger stars into a 
carbon star {Z = 0.02 standard model). The thick black 
line is the sum of all our R-merger channels, the lines with 
points are the individual channels. 

3.8. Dredge-up during the core helium flash 

We postulate that, during the core helium fiash of our 
merged objects, there must be some dredge up of carbon 
from the core into the envelope. We cannot model this 
in detail, but we can estimate the minimum amount of 
carbon, Ac, required to make C/0 = 1. We find 



Ac « Me, 



\Xo — Xc) 
1-Xc 



(1) 



where Monv is the common envelope mass and Xq and Xq 
are the abundances of carbon and oxygen initially present 
in the stellar envelope (the derivation of this formula is 
in Appendix lA.ip . In our Z — 0.02 simulations Aci2 is 
usually in the range 0.003 to 0.01 M©, as shown in figure 

m 



lOeupree ( 1986l l finds a much smaller amount of car- 
bon, around 1O~^M0, is mixed to the stellar envelope 
in two-dimensional calculations of the helium fiash. These 
simulations were necessarily of low resolution, and did not 
include rotation, which should be important in mixing ma- 
terial out of the core and th rough the hydrogen-burning 
shell ( Mengel &: Grossl . 19761 see also our section [4~4)) . 



3.9. Variation of model parameters 

The initial metallicity is the parameter which most affects 
the formation of early-R stars by the binary merger chan- 
nel. Here we consider some other parameters of our model 
which may affect the formation of early-R stars. 
— The common envelope parameter a, determines the 
fraction of the orbital energy which is used to eject 
the envelope. This most affects the initially long-period 
binary channels which go through two phases of com- 
mon envelope evolution e.g. the R3b channel. If a is 
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small the binary merges during the first common en- 
velope phase, meaning the second never happens. This 
reduces the Z — 0.02 R3 population by about 25% 
(i.e. by the fraction which are R3b stars) when chang- 
ing from a = 1 to a = 0.1. 

We have two prescriptions for A, the giant envelope 
binding energy fac tor. The first i s to u se a constant 



4. Discussion and Hope 

This raises hope of understanding these stars at last, but 
many questions remain. Here we discuss the merits of the 
different channels and consider the mixing and angular 
momentum transport which may occur during the core 
merger. 



value of 0.5, as in iHurlev et al.1 (|2Q02h - this is our 



standard mod el. Our second prescr i ption is based on 4. J. The R-chsnnels (Z = 0.02) 



stanaara moa el. Uur second prescr i ption is based 01 
the models of lOewi fc TaurisI l|200d boOlh who calcu 
lated A from detailed stellar evolution models. With 
the variable-A prescription the R-star frequency is re- 
duced by a only few per cent relative to A = 0.5. 
The amount of third dredge up in AGB stars less 
massive than 1.5 is a debated subject. Detailed 
stellar evolution models suggest no dredge up in 
these stars, but extra dredge up is required in AGB 
stars in the Magellanic clouds in order to match 



the o b served carbon sta r lum inosity function (jlben 



1981 ). Guandalini et al. ( 2006h claim we should not 



increase the dredge up, but Bonacic Marinovic et al. 
(2007, A&A submitted) confirm that extra dredge 
up is required on the basis of s-process abundances. 
Regardless, we tested the effect of increasing the 
amount of the third dredge up as a function of raetallic - 
ity according to the prescription of llzzard et al.l (|20041 ) . 
The extra dredge up increases the number of extrin- 
sic carbon stars and slightly increases the threshold 
metallicity at which the number of R3 stars is equal 
to the number of R4 stars. At solar metallicity, there 
are still no intrinsic carbon stars in our population and 
this parameter has no effect. 

Introducing the companion rei n forced attrition process 
(CRAP) of iTout fc EggletonI ()l988h has little effect 
on the number of early-R stars, reducing the num- 
ber of R3 stars by just a few per cent for large val- 
ues of the CRA P parameter (10^, as suggested in 



Tout fc EggletonI ). However, for such large values the 



R5 channel becomes important (about half as many as 
R3) because enhanced mass loss on the giant branch 
increases the number of HeWDs available for mergers. 
— Changing the initial eccentricity, e, of the binary pop- 
ulation from to 0.75, which simulates the effect of 
introducing a distribution in e, has little effect on our 
results. 



We have not considered changing the initial distributions 
such as the IMF, separation distribution etc. Such an ef- 
fort is beyond the scope of this paper. Strictly, we should 
apply a Galactic metallicity and star-formation rate distri- 
bution (perhaps as a function of height above the Galactic 
plane) rather than compare single-population, constant 
star-formation rate populations. Again, this is beyond the 
scope of this paper but we hope to examine this problem 
in the future. 



All our R channels together cannot be the equivalent of 
the observed early-R stars because we make too many by 
a factor of about ten. There are advantages and disadvan- 
tages to each channel which, in our opinion, make some 
more likely candidates than others. 

R3 The HeWD-GB mergers. The major disadvantage of 
the initially short-period R3a channel is the delay time, 
about 2.5 Myr, between the core merger and helium ig- 
nition. The equivalent in the R3b channel is less than 
1 Myr. A lengthy delay may allow transfer of material 
and/or angular momentum from the (pre-flash) core 
to the envelope, which would slow the core rotation 
and reduce the likelihood of mixing of carbon into the 
envelope during the helium flash. 
The advantage of the R3a channel is that the cores 
are of low enough mass, mostly < 0.45 M0, to be 
degenerate when helium ignites, provided they can 
lose any energy they gained from the merger pro- 
cess (which would oth erwise lift the degeneracy, as in 
Guerrero et al.l . l2004h and still remain rapidly rotat- 



ing. If Mc (the merged-core mass) is the parameter 
which determines whether the helium flash mixes ma- 
terial into the envelope then only a small range (~ 
0.02 M© e.g. 0.38- 0.40 Mq) in selects the required 
10% of R3 stars. The relatively high core mass of the 
R3b channel means its stars are many times brighter 
than the observed R stars, with logL/L0 2.5 — 3.0. 
They also span a wider effective temperature range, 
from logTcff/K = 3.55 to 3.9, than the R3a stars which 
cluster at logTcff/K ~ 3.68. Possibly these would be 
observed as J stars or as A and/or F giants. 
R3a stars, and our other R channels which involve 
signiflcant accretion from a giant companion, are en- 
hanced in ^^C and ^^N, as seen in the early-R stars 
I Dominy . 1984)- These enhancements may also be due 
to mixing across a hydrogen shell (see section 1474]) . 

R4 The binary R4 channel, equivalent to CH stars, cannot 
work at solar metallicity (see section 14.21 below) and 
does not lead to single stars, so this is definitely not 
the channel which makes the early-R stars. 

R5 Channel R5, the HeWD-HG mergers, is similar to R3 
but contributes a factor of 17 fewer stars in our stan- 
dard Z — 0.02 model. Its main advantages are that the 
number of early-R stars formed is similar to that which 
is observed and that HG stars are less compact than 
GB stars so more likely to keep the ir envelopes and 
merg e in a common envelope phase IjTaam fc Ricker . 
20061 ) . The disadvantage is that the core of a HG star 
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may not be sufficiently degenerate for helium to ignite 
soon after the merger so the core has time to lose its 
angular momentum. On the other hand, the merged 
core mass is always low (about 0.3 M©) because the 
HG star is relatively unevolved compared to an equiv- 
alent mass GB star. Assuming the merged core can 
maintain its rotation rate as the R5 star ascends the 
giant branch these may be the early-R stars. 

R6 The GB-HG mergers occur at a rate around 1/4 that 
of the R3 channel - still too many compared to the 
observations, but a better match than R3. As in R5, 
the advantage of a HG star is that it is likely to sur- 
vive the common envelope phase, but they may not 
be degenerate or massive enough to be still rotating 
when helium ignites. A fraction of R6 forms from dou- 
ble common-envelope systems and, as in channel R3b, 
these are more luminous than observed R stars. 

R7 These are the HG-HG mergers, which are very rare 
(150 times fewer than the R3 channel) so are probably 
not the early-R stars. 

R8 These AGB-helium core mergers are also rare, with 
numbers similar to R7. It is also not likely that the 
assumption of our model - that the helium core com- 
pletely mixes with the CO core of the AGB star - is 
true. In any case, these are too bright to be the early-R 
stars (but may be J/N/late-R stars). 

Rll GB-GB mergers are rare, a factor of 15 less than the 
R3 channel, but cannot be ruled out as potential early- 
R star progenitors. However, they are the binaries most 
likely to lose their envelopes when they merge. Both 
cores are degenerate and usually of low mass and often 
one member of the binary is a rejuvenated WD with 
a very thin envelope, so these are very similar to R3 
mergers. Stars in this channel span a wide range in 
luminosities (from 60 to 900 Lq) and effective tem- 
peratures (logTcff/K = 3.53 — 3.9). In a few cases the 
post-merger object has a very thin envelope, so is very 
blue (Tefi- - 40,000K). 

Rl, 2, 9 & 10 None or very few of these are made in our 
populations. 

Most likely only some fraction of our model early-R stars 
will ignite and mix carbon into the envelope, but we do 
not know which stars these are. Our population synthesis 
model predicts too many early-R stars. This is a posi- 
tive result, because we have tried to select every possible 
merger channel. It is also an interesting challenge because 
we do not know the true R-star formation channel. 



4.2. R/CH stars, metallicity and population 

Above a threshold metallicity [Fe/H] ^ —0.3 and in an ex- 
clusively old stellar population, our models suggest that 
the early-R stars should be the only carbon stars. The 
Galact ic bulge, which is o ld and of (super-)solar metal- 



have been identified as early-R stars in the surveys of 
Azzopardi. Lequeux fc RebeirotI (|1985t ): . This is in con- 



licity ( Zoccali et al . 2003f ). is an obvious target for R 
star surveys. In fact, large numbers of low-luminosity 
carbon stars have already been found in the bulge and 



trast to th e very few o r even total lack of N-type stars in 
the bulge (jRichl . I 119981 ). in agreement with our results. 

At metallicities typical of the Galactic halo, the num- 
ber of CH stars vastly exceeds the number of early-R stars. 
That is not to say there are no early-R stars in the halo, 
because according to our models there should be. With 
regard to the detection of early-R stars in the halo, there 
are several reasons it may be difficult: 

1. There should be many more CH stars than early- 
R stars in the halo. At [Fe/H] = —2.3 we predict 
R3/R4 ^ 7% which means the early-R/CH ratio will 
be much smaller (R4 only includes core helium burn- 
ing CH stars, the true CH population will contain both 
these and giants and possibly dwarfs) especially if we 
reduce our figure by a factor of ten to agree with the 
Hipparcos observations. 

Halo stars are rare compared to disk stars, so there are 
simply not as many population-II early-R stars. 
CH stars are likely to be selected by colour, which may 
omit early-R stars. 

Giant branch CH stars are somewhat brighter than 
CHeB stars, so any survey aimed at CH stars (rather 
than C dwarfs) must be quite sensitive to relatively 
dim stars in order to detect an early-R star. 
Is it possible to distinguish easily between early-R 
and (core-heliura burning) CH stars? According to 
Abia et al. 1 20031 ) the only differences are the presence 



2. 



3. 



5. 



of s-process isotopes and binarity, data which are often 
difficult to obtain. As such, early-R stars already found 
in the halo could have been misclassified as CH stars. 
A recent survey has detected a few halo early-R stars 
and gives a me thod for d i fferen tiating between early-R 
and CH stars I Goswariiil . 2005f l. 

The IWesterlund et al.l (11995!) survey of the Large 



Magellanic Cloud (see also iRebeirot et al.l . 119931). which 



has a sub-solar metallicity [Fe/H] ~ —0.3 (jCioni et al 



2006al ). finds a menagerie of N, R and J type stars, over a 



wide range of luminosities, in qualitative agreement with 
our model population. We should caution, however, that 
the LMC probably contains a s ignificantly younger popu- 
lation of intrinsic N type stars i Cioni et al. . 2006b[ ) which 
makes a statistical comparison difficult because our pop- 
ulation model does not include young stars. 

We also have not considered the impact of the 
latest determinations of the solar oxyge n abun- 



dance tAllende Prieto. Lambert fc Asplundl . 
Asplund. Grevesse fc Sauval 2005l l which decreases 



2001 



the amo unt of oxygen b y a factor of two relative to 
[Anders fc Grevesse 1 19891 ). This reduces our estimate of 
the mass of carbon which must be mixed into the envelope 
during helium ignition, Aci2, and shifts the metallicity 
at which CH stars form, but does not significantly change 
our conclusions. 

Finally, we are not sure why Galactic early-R star s 
are located only in the thick disk i Bergeat et al. . 2002bl ). 
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Some of their HCl-3 (R-type) stars are of low radial ve- 
locity, so could be members of the thin disk. Alternatively, 
the initial binary fraction, and hence number of mergers, 
may be higher in the th i ck dis k compared to the thin disk 



I Grether fc Lineweaver . 20061 ). 



4.3. Progenitors observed 

The progenitors of our early-R stars are binaries in which 
one star is a WD, the other a blue straggler - one such 
star has been observed (jLandsman et al.1 . Il997l ) . Recently, 
models of MS+WD blue st raggler systems were con- 
structed by lTian et al. ( 2006h but they cannot take their 
calculations through the final common envelope phase. 
Regarding the MS+WD phase, however, they find quite 
similar results to ours despite their use of full stellar evo- 
lution rather than synthetic code. 

Blue stragglers (BSs) are most easily identified in glob- 
ular clusters (GCs) because the stars in a particular cluster 
are of the same age so BSs are easily seen to be brighter 
and bluer than the tip of the main sequence. However, 
binary mergers may be rare b ecause of the low b inary 
fraction in the core of clusters ( Ivanova et al 112005) - or 
perhaps the binary fraction is low because of mergers. To 
search for progenitors of our R-star mergers it would be 
best to search for HeWD-BS binaries in high metallicity, 
old GCs, such as those in the Galactic bulge e.g . NGC 



6553 or 6528, (jBarbuv et al.1 . 12004 : IZoccaH et al.1 . 120041 ) 



with metallic i ties [ Fe / H] 
i Llovd Evanl . ll986h . 



-0.2, or the LMC clusters 



The FK C om stars (|Bopp fc Rucinski 
Bopp fc Ste^ . Il98lh may be post-common 
lope mergers: they are single, ra pidly rotating 
to or even beyon d brea k-up e.g. iHuenemoerder et al 



enve- 
(close 



1993t lAvres et al.1 l2006l ). G/K-type giants which are 
thought to have e volve d through a binary merger 



i Weltv fc Ramsevl. 



10- 



.x.ov.... ,19941). Th ev a re rare, about 

pc^'^ I Collier Cameronl . 1982; Hagen fc Stencej 



198,'tI ). which is quite similar to the R-star space density. 



4.4. Merger mixing and modelling 

In our merger scenario aX Z = 0.02 the amount of carbon 
which is required to be mixed from the core into the enve- 
lope during the helium fiash is 0.003 — 0.01 Mq , assuming 
C/0 = 1 in the R-star envelope. The analysis of iDominv 
1 19841 ) shows carbon and nitrogen are similarly enhanced 
in early-R stars, with the ^^C/^^C ratio just above the CN- 
cycle equilibrium value of 4. These observations strongly 
suggest there is some hydrogen burning of the carbon-rich 
material from the helium fiash en route to the surface of 
the early-R star. This would mean we have underestimated 
the amount of carbon dredge up by a factor of about two 
because half the carbon is converted to nitrogen. However, 
because some ^^C and ^"'N may come from the deposition 
of material stripped from the primary star during the first 



period of mass transfer, for some channels (e.g. R3a) we 
require less CN burning. 

As discussed above, there may be a delay between the 
merger and helium ignition. If the core can settle to hy- 
drostatic equihbrium it should be possible to model it 
with a one-dimensional stellar evolution code which in- 
cludes rotation and angular momentum tran sport (e.g. 
Langer , 1998t Heger. Langer fc Wo oslevl. 2000l and the se- 
ries of papers by Meynet fc ^Maedej, Il997 . onward). The 
SPH models of FCuerrero et al.l |200J) suggest that merged 
cores settle down rapidly (minutes) compared to their nu- 
clear burning timescale (Myr) so the problem is tractable. 
Certainly, some one-dimensional exploration of the prob- 
lem would be useful. 

This leads us to the realm of two and three di- 
mensional modelling of the core merger and associated 
mixing. Two dir n ensional simulations without rotation 
(ICole fc Deupreel . Il980l . Il98ll : lOeupree fc Wallace! . Il987l : 
Deupreel . ll996f ) do not show significant mixing of material 
from the core to the surface - yet we have good grounds, 
based on our population model and a lack of alternatives, 
to suggest that it does happen. The only way forward may 
be a full three-dimensional simulation of a merger of two 
HeWDs inside a hydrogen-rich giant envelope including 
nuclear burning. This is some way off, and might not even 
be possible with current methods such as SPH, but is per- 
haps not impossible in the near future using codes such as 
Djehuty l|Dearborn et all . l2Q06t ). 



198 It 4.5. Angular momentum 



In appendix IA.2I we show that when two HeWD cores 
merge the velocity of the resultant object exceeds its 
brea kup velocity by a factor of abo ut two, in agreement 
with GourgouHatos fc Jeffery 1 20061 ). There must be some 
form of angular momentum transport out of the merged 
core, perhaps by magnetic fields coupling the core to the 
envelope or shear mixing across the core -envelope bound- 
ary ( it has recently been suggested by ISoker fc Harpad . 
20071 that shear mixing may mix helium-rich, but not 
carbon-rich, material out of the core of a red giant). 

We show in appendix IA.3I that if we conserve angular 
momentum and mass for the whole star, after the merger 
it is rotating at about its breakup velocity. However, it 
is not the case if some mass is lost during the common 
envelope phase as this may remove enough angular mo- 
mentum to stabilise the star ( Davies. Benz fc HillfJ . ll99ll 



iLovekin fc Sillsl . l2003f ). It is still rotating at a significant 
fraction of its breakup velocity, ofte n many tens of k m . 
This contradicts the observations o f lMcClurd l|l997l ) which 
show that early-R stars are not rapidly rotating. There 
must be some angular momentum loss, such as magnetic 
braking, which slows the early-R star, especially because 
it contracts after helium ignites at the tip of the giant 
branch. 
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4.6. After the R-star. . . 

The early-R stars are CHeB stars but eventually they 
will exhaust their supply of heHum and ascend the gi- 
ant branch again as ^^C-, ^^C- and ^^N-rich AGB stars. 
They may be observed as J-type stars and indeed, per- 
haps early-R stars are the proge nitors of the J stars, as 
suggested by iLloyd Evans (|l986h . If this is the case then 
J stars should also be single stars, but it is not clear if 
this is tru e. J stars are appare ntly 15% of carbon-rich gi- 
ants stars ( Abia &: Isern . 200Cll l and further work regarding 
their formation rate will determine if early-R stars are the 
progenitors of some of them. 



5. Conclusions 

We have shown by means of a binary population synthesis 
technique that binary mergers are a viable channel for the 
formation of the early-R stars. Our models naturally re- 
produce the properties of the early-R stars, namely that 
they are single, K-type stars which show enhanced ^^C, 
^^C and ^^N without s-process or iron enhancement rela- 
tive to the sun. 

We estimated the number of early-R stars from differ- 
ent binary merger channels involving objects with helium 
cores. The most common merger is that of a helium white 
dwarf with a red giant star and this makes about ten times 
as many stars as we require if we are to match the early-R 
to red clump number ratio to the observed value. We did, 
however, assume that the helium core of every merger is 
rotating and mixes carbon into the envelope upon helium 
ignition. It is more likely that a small fraction of our stars 
ignite heHum while rotating rapidly enough to mix carbon 
into their envelope. 

Other possible channels include mergers of a red giant 
with a Hertzsprung gap star, a heHum white dwarf with 
a Hertzsprung gap star or the merger of two red giants. 
There must be other parameters which control whether a 
merged star becomes an early-R star, such as the core mass 
or core-rotation rate. Our model early-R stars are rapidly 
rotating core-heHum burning giants so must undergo some 
kind of rapid magnetic braking to slow them down, be- 
cause observed early- R stars are not rota ting quick ly. 

As pointed out bv lPominvl l|l984[ ) and lMcClurl l|l997l l 
it is difficult to envisage an alternative evolutionary path- 
way which leads to the early-R stars. Ten years on, we fi- 
nally have a potential explanation for their existence and 
their number. It is our hope that this will stimulate fur- 
ther research in the area of helium-core mergers inside 
common envelopes because for once we think we under- 
stand the outcome of the merger process. 



Acknowledgements. Our work on the early-R stars com- 
menced at the "Nucleosynthesis in Binary Stars" work- 
shop at the Lorentz centre, Leiden, 2005. RGI is sup- 
ported by the Nederlands Organisatie voor Wetenschappelijk 
Onderzoek (NWO). He thanks Armagh Observatory and 
Monash University for funding extended visits. He is grate- 
ful to the numerous people who have shared their ideas with 



him over the past year on the subject of the paper, es- 
pecially Evert Glebbeek, Onno Pols and Pilar Gil-Pons for 
reading the manuscript. CSJ acknowledges financial support 
from the Northern Ireland Dept. of Culture Arts and Leisure 
(DCALNI), the UK Particle Physics and Astronomy Research 
Council (PPARC) and the Netherlands Research School for 
Astronomy (NOVA). Other aspects of this research were sup- 
ported by the DEST and the Australian Research Council. This 
research has made use of the SIMBAD database, operated at 
CDS, Strasbourg, France. 

References 

Abia C, Dominguez I., GaHino R., Busso M., Masera S., 
Straniero 0., de Laverny P., Plez B., Isern J., 2002, ApJ, 
579, 817 

Abia C, Dommguez I., Gallino R., Busso M., Straniero 

0. , de Laverny P., WaHerstein G., 2003, PubHcations of 
the Astronomical Society of AustraHa, 20, 314 

Abia C., Isern J., 2000, ApJ, 536, 438 
Allende Prieto C, Lambert D., Asplund M., 2001, ApJ, 
556, L63 

Anders E., Grevesse N., 1989, Geochimica et 
Cosmochimica Acta, 53, 197 

Asplund M., Grevesse N., Sauval A. J., 2005, in Barnes 
III T. G., Bash F. N., eds, ASP Conf. Ser. 336: 
Cosmic Abundances as Records of Stellar Evolution and 
Nucleosynthesis The Solar Chemical Composition, p. 25 

Ayres T. R., Harper G. M., Brown A., Korhonen H., Ilyin 

1. v., Redfield S., Wood B. E., 2006, ApJ, 644, 464 
Azzopardi M., Lequeux J., Rebeirot E., 1985, A&A, 145, 

L4 

Barbaro G., Dallaporta N., 1974, A&A, 33, 21 
Barbuy B., Melendez J., Ortolani S., ZoccaH M., Bica E., 
Renzini A., HiH V., Momany Y., Minniti D., Rich M., 
2004, Memorie deHa Societa Astronomica ItaHana, 75, 
398 

Barnbaum C, Stone R. P. S., Keenan P. C, 1996, ApJS, 
105, 419 

Baumert J. H., 1974, ApJ, 190, 85 

Bensby T., Feltzing S., Lundstrom I., 2004, A&A, 421, 
969 

Bergeat J., Knapik A., Rutily B., 2002a, A&A, 390, 967 
Bergeat J., Knapik A., Rutily B., 2002b, A&A, 385, 94 
BesseH M. S., Wood P. R., Lloyd Evans T., 1983, MNRAS, 
202, 59 

Blanco V. M., 1965, in Blaauw A., Schmidt M., eds. 
Galactic Structure Distribution and Motions of Late- 
Type Giants, p. 241 
Bopp B. W., Rucinski S. M., 1981, in Sugimoto D., 
Lamb D. Q., Schramm D. N., eds, lAU Symp. 
93: Fundamental Problems in the Theory of Stellar 
Evolution The Rapidly Rotating Giants of the Fk- 
Comae Type. p. 177 
Bopp B. W., Stencel R. E., 1981, ApJ, 247, L131 
Cannon A. J., Pickering E. C, 1918, Annals of Harvard 

CoHege Observatory, 91, 1 
Cannon R. D., 1970, MNRAS, 150, 111 
Cannon R. D., Lloyd C, 1969, MNRAS, 144, 449 



R.G. Izzard et al.: Origin of the early-type R stars: a binary-merger solution to a century-old problem 



15 



Cioni M.-R. L., Girardi L., Marigo P., Habing H. J., 2006a, 
A&A, 448, 77 

Cioni M.-R. L., Girardi L., Marigo P., Habing H. J., 2006b, 

A&A, 452, 195 
Cole P. W., Deupree R. G., 1980, ApJ, 239, 284 
Cole P. W., Deupree R. G., 1981, ApJ, 247, 607 
Collier Cameron A., 1982, MNRAS, 200, 489 
Dahn C. C, 1964, PASP, 76, 403 

Davies M. B., Benz W., Hills J. G., 1991, ApJ, 381, 449 
Dearborn D. S. P., Lattanzio J. C, Eggleton P. P., 2006, 

ApJ, 639, 405 
Deupree R. G., 1986, ApJ, 303, 649 
Deupree R. G., 1996, ApJ, 471, 377 
Deupree R. G., Wallace R. K., 1987, ApJ, 317, 724 
Dewi J. D. M., Tauris T. M., 2000, A&A, 360, 1043 
Dewi J. D. M., Tauris T. M., 2001, in Podsiadlowski P., 

Rappaport S., King A. R., D'Antona F., Burderi L., eds, 

ASP Conf. Ser. 229 p. 255 
Dominy J. F., 1984, ApJS, 55, 27 
Dray L. M., Tout C. A., 2003, MNRAS, 341, 299 
Eggen 0. J., 1972, MNRAS, 159, 403 
Eggleton P. P., 1983, ApJ, 268, 368 

Eggleton P. P., Kiseleva-Eggleton L., 2002, ApJ, 575, 461 
Faulkner D. J., Cannon R. D., 1973, ApJ, 180, 435 
Fleming W., Pickering E. C, 1896, Harvard College 

Observatory Circular, 9, 1 
Fleming W., Pickering E. C, 1908, Harvard College 

Observatory Circular, 145, 1 
Goswami A., 2005, MNRAS, 359, 531 
Gourgouliatos K. N., Jeffery C. S., 2006, MNRAS, 371, 

1381 

Grether D., Lineweaver C. H., 2006, ArXiv Astrophysics 

e-prints, astro-ph/0612172 
Guandalini R., Busso M., Ciprini S., Silvestro G., Persi 

P., 2006, A&A, 445, 1069 
Guerrero J., Garcia-Berro E., Isern J., 2004, A&A, 413, 

257 

Hagen W., Stencel R. E., 1985, AJ, 90, 120 
Harm R., Schwarzschild M., 1966, ApJ, 145, 496 
Heger A., Langer N., Woosley S. E., 2000, ApJ, 528, 368 
Hollowell D., Iben I. J., Fujimoto M. Y., 1990, ApJ, 351, 
245 

Huenemoerder D. P., Ramsey L. W., Buzasi D. L., Nations 

H. L., 1993, ApJ, 404, 316 
Hurley J. R., Pols O. R., Tout C. A., 2000, MNRAS, 315, 

543 

Hurley J. R., Tout C. A., Pols 0. R., 2002, MNRAS, 329, 
897 

Iben I., Renzini A., 1983, ARA&A, 21, 271 
Iben Jr. I., 1981, ApJ, 246, 278 
Iben I. J., 1990, ApJ, 353, 215 
Ishida K., 1960, PASJ, 12, 214 

Ivanova N., Belczynski K., Fregeau J. M., Rasio F. A., 

2005, MNRAS, 358, 572 
Izzard R. G., Dray L. M., Karakas A. I., Lugaro M., Tout 

C. A., 2006, A&A, 460, 565 
Izzard R. G., Tout C. A., Karakas A. I., Pols 0. R., 2004, 

MNRAS, 350, 407 



Karakas A. I., Lattanzio J. C, Pols O. R., 2002, PASA, 
19, 515 

Keenan P. C, Heck A., 1994, Revista Mexicana de 

Astronomia y Astrofisica, vol. 29, 29, 103 
Keenan P. C, Morgan W. W., 1941, ApJ, 94, 501 
Kippenhahn R., Kohl K., Weigert A., 1967, Zeitschrift fur 

Astrophysik, 66, 58 
Knapp G., Pourbaix D., Jorissen A., 2001, A&A, 371, 222 
Kroupa P., Tout C, Gilmore G., 1993, MNRAS, 262, 545 
Landsman W., Aparicio J., Bergeron P., Di Stefano R., 

Stecher T. P., 1997, ApJ, 481, L93 
Langer N., 1998, A&A, 329, 551 

Lattanzio J., Dearborn D., Eggleton P., Dossa D., 2006, 

ArXiv Astrophysics e-prints, astro-ph/0612147 
Livio M., Soker N., 1988, ApJ, 329, 764 
Lloyd Evans T., 1986, MNRAS, 220, 723 
Lodders K., Fegley Jr. B., 1998, Meteoritics and Planetary 

Science, 33, 871 
Lovekin C, Sills A., 2003, Canadian Undergraduate 

Physics Journal, II, 15 
McClure R. D., 1997, PASP, 109, 256 
Mengel J. G., Gross P. G., 1976, Ap&SS, 41, 407 
Meynet G., Maeder A., 1997, A&A, 321, 465 
Paczynski B., Tremaine S. D., 1977, ApJ, 216, 57 
Rebeirot E., Azzopardi M., Westerlund B. E., 1993, 

A&AS, 97, 603 
Regos E., Tout C. A., 1995, MNRAS, 273, 146 
Rich R. M., 1998, in Sofue Y., ed., lAU Symp. 184: 

The Central Regions of the Galaxy and Galaxies The 

Galactic bulge (Review), p. 11 
Rufus W. C, 1916, Publications of Michigan Observatory, 

2, 103 

Sanford R. F., 1944, ApJ, 99, 145 

Scalo J. M., 1976, ApJ, 206, 474 

Secchi A., 1868, MNRAS, 28, 196 

Shane C. D., 1928, Lick Observatory Bulletin, 13, 123 

Soker N., Harpaz A., 2007, ArXiv Astrophysics e-prints. 



astro-ph/0701528 



Stancliffe R. J., Glebbeek E., Izzard R. G., Pols O. R., 

2007, A&A, 464, L57 
Stephenson C. B., 1973, Publications of the Warner & 

Swasey Observatory, p. 1 
Taam R. E., Ricker P. M., 2006, ArXiv Astrophysics e- 



prints, |astro-ph/0611043 ) 
Tian B., Deng L., Han Z., Zhang X. B., 2006, A&A, 455, 
247 

Tout C. A., Aarseth S. J., Pols 0. R., Eggleton P. P., 1997, 

MNRAS, 291, 732 
Tout C. A., Eggleton P. P., 1988, MNRAS, 231, 823 
Vandervort G. L., 1958, AJ, 63, 477 
Wallerstein G., Knapp G. R., 1998, ARA&A, 36, 369 
Welty A. D., Ramsey L. W., 1994, ApJ, 435, 848 
Westerlund B. E., Azzopardi M., Breysacher J., Rebeirot 

E., 1995, A&A, 303, 107 
Zamora O., Abia C, Plez B., Dommguez I., 2006, Memorie 

della Societa Astronomica ItaHana, 77, 973 
Zoccali M., Barbuy B., Hill V., Ortolani S., Renzini A., 

Bica E., Momany Y., Pasquini L., Minniti D., Rich 



16 



R.G. Izzard et al.: Origin of the early-type R stars: a binary-merger solution to a century-old problem 



R. M., 2004, A&A, 423, 507 
Zoccali M., Renzini A., Ortolani S., Greggio L., Saviane 
I., Cassisi S., Rejkuba M., Barbuy B., Rich R. M., Bica 
E., 2003, A&A, 399, 931 

Appendix A: Appendix 

A.l. Carbon Dredge Up 

We define abundances by mass fraction to be X prior to 
the merger, Y post-merger. We then require 



41c 
31o 



> 1 



(A.l) 



to be a carbon star. The common envelope is of mass M 
and the amount of carbon required from the core is Ac 
(a similar mass of material is mixed back down into the 
core). We neglect the minor isotopes ^^C, ^''O and ^^O 
because their abundances are small (at most 1/4 of C is 
^■^C). The final mass of carbon in the envelope is 



FcMenv = Xc(M- Ac) + Ac 
which gives 



Xc(Af- Ac) + Ac 
M 



(A.2) 



(A.3) 



If we assume the oxygen abundance is unchanged, i.e. 
Fo = -'^O, and there is no CN cycHng of dredged up car- 
bon, we have 



(A.4) 



which with the previous expression and some algebraic 
manipulation gives us 



Ac > 



|Xo — Xc 
1-Xc 



M. 



(A.5) 



If a fraction / of the dredged-up carbon is converted to ni- 
trogen as it passes through the hydrogen shell, the amount 
of dredged-up carbon and nitrogen is 



A 



CN 



> 



1 l-'^o 



Xr 



-M. 



(A.6) 



A.2. Merged Core Angular Momentum 



Here we consider the angular momentum in the early-R 
star core as a result of the merger. When the least massive 
helium core, of mass mc2, enters RLOF, the orbital period 
is given by Kepler's law 



P 




TTlcl -I- TOc2 



(A.7) 



where ol is the separation and TOcI is the mass of the more 
massive core. The equivalent orbital angular velocity is 



then = 2tt/P. We can estimate the angular momentum 
in the twin-core system J by neglecting the spins of the 
stars such that 



J — Jcl + Jc2 + Jovh,! 



(A., 



where Jci,2 are the spin angular momenta of the stars and 
H is the reduced mass. Then we have, by substitutions of 
the expression for ft, 



Jo 



rrici + mc2 



rial 



TOcimc2 



mci -I- mc2 



(A.9) 



Now, if we conservatively assume that no angular momen- 
tum is lost from the cores during the merger then we can 
calculate the orbital angular velocity lo of the merger prod- 
uct, 



Jo: 



knir 



TOcimc2\/G'aL 
/cr2 (mci + TOc2)' 



(A.IO) 



where rric = md -I- mc2 and Tc are the mass and radius of 
the merged core respectively and k is its radius of gyration. 
The velocity at the surface of the merged core is 



mcimc2VGa^ 



kvc {rrici + mc2)^ 
We can easily compare this to the break-up velocity 



(A.ll) 



to see 



Grric 



TOcimc2 



mc2) 



(A.12) 



(A.13) 



The separation at 
max frn/ffq), Tr-?./ fiH q)] 



RLOF is OL = 

_. where / is given by the 

function of lEggleton 1 1983f ). and given md and TOc2 we 
obtain the core radii from the formulae of iTout et aL 
()l997l ). We set k = 0.21 as for an n = 3/2 polytrope. For 
0.1 < TOd = TOc2 < 0.5 the ratio 1.9 < vy/vb ^ 2.7 so 
there must be some form of outward angular momentum 
transport in order for the cores to merge. 

A.3. Rotational Velocity of the Post-CE Star 

We can derive the maximum velocity of the star which re- 
sults from common envelope evolution by conserving an- 
gular momentum during the merger. At the beginning of 
the CE phase the binary system has a total angular mo- 
mentum J = Jorb + Ji + J2 where Ji and J2 are the spin 
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angular momenta of the stars and Jorb is the orbital angu- 
lar momentum. The final angular velocity lo is calculated 
from a solid-body approximation, assuming angular mo- 
mentum is conserved during the merger, 

J = (fc2Monv + hAQR'^uj (A.14) 

where Mcnv , Mc, R are the envelope and core mass and 
the radius of the merged star respectively. The constant 
fes = 0.21 as previously, but ^2 depends on the struc- 
ture of the giant star (it is about 0.15 — 0.2). We then 
solve for the velocity v = ojR and compare to the breakup 
velocity vb (c.f. section [A.2p . For a typical R3a system, 
with initial parameters Mi = 1.25 Mq, M2 = 0.5 Mg, 
P = 1.6 days, Z = 0.02, acE = 1, Aqe = 0.5, e = 0, no 
extra CRAP, we find the merged stellar mass is 1.34 Mq 
(0.05 M0 is lost during the common envelope phase, we 
have not taken into account the loss of angular momen- 
tum due to this mass) with a core of mass 0.41 Mq and 
k2 = 0.143, hence v = 125km/s, about 95% of the break 
up velocity (132km/s). The ratio reduces to 60% if we 
assume the angular momentum in the cores does not con- 
tribute to the rotation of the envelope, which is required 
for the core to be rapidly spinning when it ignites helium. 
If we attempt to take into account the angular momentum 
that is lost in the 0.05 Mq of material ejected during the 
common envelope phase, presumably at approximately the 
breakup velocity so A J — vbR x 0.05 Mq, the ratio v/vb 
drops further to 28% {v = 37km/s). It would seem these 
systems are, in general, quite rapidly rotating (at least 
tens of km/s), but are not all at the breakup velocity. 



